ACE Sept. 19'i2 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



WARTIME REPORT 

ORIGINALLY ISSUED 

September 19^2 as 
Advance ConTidential Report 

MAEfSIS OF HEAT AND COMPRESSIBILITY EFFECTS IK IREEBHAl 
FLOW SSTSDBB A1JD HIGH-SPEED TESTS 
OF A RAM-JET SYSTEM 
By John V. Becker and Donald D. Baals 

Langley Memorial Aeronautical Laboratoiy 
Langley Field, Va. 



NACA 



WASHINGTON 



NACA WARTIME REPORTS are reprints of papers originally issued to provide rapid distribution of 
advance research results to an authorized group requiring them for the war effort. They were pre- 
viously held under a security status but are now unclassified. Some of these reports were not tech- 
nically edited. All have been reproduced without change in order to expedite general distribution. 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



ADVANCE CONFIDENTIAL REPORT 



ANALYSIS OF HEAT AND COMPRESSIBILITY EFFECTS IN INTERNAL 
FLOW SYSTEMS AND HIGH-SPEED TESTS 
OF A RAM-JET SYSTEM ' 
By JOHN V. BECKER AND DONALD D • BAALS 



SUMMARY 

An analysis has been made "by the NACA of the effects of heat and 
compressibility in the flow through the internal systems of aircraft. 
Equations and charts are developed whereby the flow characteristics at 
key stations in a typical internal system may "be readily obtained. 
The effects of compressibility (density change) were fully accounted 
for and were found to be important at present-day flight conditions, 
particularly at high altitudes. Particular attention was given to 
the .flow changes across radiators and air-cooled engines-. It is 
shown -that very high Mach numbers will be attained within the fins of 
existing air-cooled engines if required cooling is to be maintained at 
.. altitude . The density decrease and the pressure drop across radiators 
: and, engines are shown to be considerably greater than the values 
computed on the basis of an incompressible flow.' The equations show 
that when a Mach number of 1.0 is reached at any station of the internal 
system, the flow becomes choked and any further reduction of outlet 
pressure will not increase the flow rate . The addition of heat had a 
throttling effect on the flow rate and led to choking at lower rates 
of flow. . 

Tests of a 13-inch diameter ram- jet type of propulsion system 
using an electrical heating device to. add heat to the internal flow 
have been made in the Langley 8-foot high-speed wind/ tunnel throughout 
a. Mach number range of 0.2 to 0.75- The results showed that the 
measured thrusts, and thermal efficiencies closely approached the ideal 
values for this type of propulsion system. Measurements of the flow 
details at various stations of the internal system were compared with 
calculations based on one -dimensional compressible flow theory, and 
excellent agreement, was obtained. 

In regard to the propulsive effect' due to heating the internal- 
flow, the thermal efficiency of the energ7 -recovery cycle was show 
to depend mainly on the* ratio of the static pressure In the duct 
behind the radiator to the stream static pressure. To a lesser extent, 
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the thermal efficiency was governed by the Mach number of the duct flow 
and "by the rate of heating; increases in either of these quantities 
caused decreases in the thermal efficiency . 



INTRODUCTION 

I 

Tlie operating speeds and the altitudes attained by modern aircraft ^ 
have advanced to the point where the use of simple incompressible-flow 
relations in the design of internal ducts and cooling systems results 
in large errors in the- determination of the internal-flow characteristics. 
A more comprehensive 'analysis esrploying compressible -flow relations , 
therefore, is required in order to compute accurately the conditions 
at each station within the duct, the pressure drop across the radiator, 
■the internal mass flow, and. the variation of the internal-flow charac- 
teristics with the addition of heat. The effects of the addition of 
heat to the internal flow are also of importance because of the 
possibility of obtaining appreciable thrust power from the internal 
system. 

References 1 and 2, which are representative of the numerous 
theoretical analyses comprising most of the previous work in this 
field, were concerned mainly with evaluating the net force resulting 
from the flow through the engine-cooling system. In these earlier - 
analyses it was usually not necessary to compute the internal-flow 
characteristics ' in great detail and it was possible to neglect, in part, 
the effects of compressibility (density changes). A more comprehensive 
general discussion of the cooling- system. air flow was given in 
reference 3 "but the details of the flow relations involved were not 
shown . : 

In the present investigation a theoretical analysis of the flow 
details at key stations in a typical system was made. The effects of 
compressibility and of heating were accounted for up to Mach numbers of 1.0. 
Results of this analysis are. presented in part I of the investigation 
in the form of charts and working equations- that can be readily applied 
in practice* 

Of great interest are the designs of internal flow systems that 
employ much higher rates of heating than are normally found in the 
cooling systems alone. The effectiveness of such designs depends on 
the thermal efficiency of the process whereby heat is converted into 
thrust power. For this reason, one of the' -principal purposes of this . 
investigation was t9 study analytically the energy -recovery process 
and to determine by experiment how closely the ideal thermal efficiency 
could be approached in a typical system. 
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An experimental investigation of an idealized ram- jet propulsion 
system was conducted in the 8-foot high-speed wind tunnel at the 
Langley Memorial Aeronautical Laboratory in March 19^1. The object 
of the present tests was to provide experimental data to establish 
the validity of the analytical results of part I, particularly with 
reference to the conversion of heat energy. to thrust power. The results 
of these tests are presented in part II. 

The body used in this investigation was designed to confomi closely, 
with the duct arrangement assumed in. the analysis given in part I. 
Heat was added by a specially constructed radiator causing "negligible 
blockage of the duct flow and insignificant friction loss. This manner 
of adding heat minimized the fire hazard in the wind tunnel and 
eliminated the combustion problem that would have existed if fuel. were 
burned in the -duct* The amount of heat added (l6o kw max.) produced 
a thrust equal to about kO percent of .the drag of the test body . 
at M = 0.75. This rate of heating , although lower than the rates 
obtainable by combustion of fuel/ was. large enough to permit accurate 
me.as-ureaaents- of. ..the. ..thrust. ..and .other effects .of . heating to be obtained 
and thus to provide the basis .for significant comparisons with the 
analytical investigation. 



SYMBOLS 



area, square, feet 



a 



velocity of sound, feet per second 




internal-drag coefficient 




AC-p internal -drag- coefficient increment 




drag coefficient .of radiator tube 




r 




specific heat at constant pressure (for air 0.21* Btu/lb/°F) 
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D net drag due to internal flow, pounds. 

AD : internal-drag increment, pounds 

D f . . drag force due to friction within the radiator 
tubes , pounds * 

D r total drag force acting on the radiator, pounds 
d *• duct diameter t - inches 

F maximum fuselage cross-section area, square feet 

(1.009 for model tested) 

g , acceleration of gravity (32.3 ft/sec 3 ) 

H heat added, Btu T>er second 

h , total pressure, pounds per square foot absolute 

Ah total-pressure loss, pounds per square foot 

J mechanical equivalent of heat (778 ft-lla/Btu) 

K.E. kinetic energy, foot-pounds 

L. over-all length of fuselage, inches 

I distance from nose, inches 

M Mach number (v/a) 

m mass-flow rate, slugs per second 

/ _ \ . 

P pressure coefficient I- — 

p static pressure, pounds per square foot absolute 

Ap static-pressure, decrease , pounds per square foot 

Q, quantity of flow, cubic feet per second , > 

q dynamic pressure, pounds per square foot X pVT . 

R gas constant, feet per °F (for air, 53.3) 



R radius , inches 

T temperature, °F absolute 

I thrust, pounds 

o — 

At temperature change, F 

V 0 free-stream velocity, feet per second 

v velocity within duct, feet per second 

x distance from leading edge of respective sections 

inches 

Y ratio of specific heats (for air, 1.40) 

€ fi heat-cycle efficiency (equation (31)) 

e M . mechanical-efficiency factor (equation (33)) 

€j propulsive efficiency (c^Cjj) (equation (32)) 

p. density, slugs per cuMc foot 

n d dif fuser- ef f ici ency 

(h - p M s M 4 M 6 

1+ri compressibility factor ^— - — = 1 + — + ^goo 

Subscripts : 

0 to 5 stations in internal flow system shown in figure 
c test condition without heat 

d value at any station in duct 

f friction component 

1 low speed, incompressible-flow condition 

r condition . across the radiator or heating device 

r 2 within the radiator at the tube entrances 

r^ within the radiator at the tube exits 

t condition within the radiator tubes 
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A prime after a symbol indicates the condition with- 
out heat "but with same mass flow as with heat. 

I, ANALYSIS OF INTERNAL FLOW SYSTEMS 
Calculation of the Flow Characteristics at Key 
Stations in a Typical Internal Flow System 

In the design of efficient aircraft, the flow char- 
acteristics (static pressure, density, and velocity) must 
he computed at several key. stations in each internal flow 
system. The ohjects of such calculations are to deter- 
mine the state of the air entering the radiators, the 
pressures available for cooling,, the flow changes across 
the heating device, the net drag or propulsive force due 
to the internal flow, and the required sizes of inlet and 
outlet openings. It will "be shown that under present-day - 
operating conditions compressibility effects and the sec- 
ondary effects of heating cannot be neglected in the cal- 
culations. 

Methods will he presented for analyzing the internal 
flow system that will account for the effects of heating 
and compressibility and will yet be simple to use. Elim- 
ination of much of the mathematical work will be effected 
through solution of the more cumbersome equations in chart 
form. The methods presented are generally applicable to 
any internal flow system, either with or without the addi- 
tion of heat or mechanical energy (internal blower), and 
for any rate of internal flow. 

One-dimensional flow and a uniform velocity distribu- 
tion across the duct, is assumed throughout the analysis*,, 

The theory of. the process whereby the heat energy 
added to the internal flow is converted into mechanical 
energy is described under a subsequent subheading of part 
I. 

Initial conditions . - The initial quantities assumed 
to be known by the designer, are as. follows: 

(l) The rate ■ at which heat is to be dissipated H , 
Btu per second 
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(2) Th<e pressure drop for the cold radiator (Ap/qg^ 

obtained from low-speed tests 1 

(3) The mass-flow rate required by the system m, 

slugs per second 

(4) The speed of flight ? 0 and the atmospheric 

conditions p Q , p Q , and .T 

(5) The duct areas available at all stations in the 

internal system 

Station 0,- In the method to be described, the inter- 
nal flow is considered as a stream tube starting far ahead 
of the airplane Where the velocity in the tube (relative 
to the airplane) is V 0 and the air conditions are those 
of the free stream (fig. l). The. area of the tube at 
this station is determined by the required rate of inter- 
nal flow; that is, 



KO V 0 

The area A 0 is noted to be a measure of the mass 
flow and bears a simple relation to the flow coefficients, 
Q, g 

FV 0 a P 0 F ^7* If ^ equals the volume rate of flow 
at any station d, then 

0. _ m_ _1_ Po fo 
The mass-flow coefficient C is 



C - m - — 

" Po FV o ' 

For isentropic flow, the pressure, the density, and the 
temperature at any station in the duct for any rate of 
internal flow can be related to the known initial con- 
ditions through the area ratio A 0 /A d and the Mach number 
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U Q zz -2- f where a Q = 49.0 sjrj. Charts for obtaining 
a o 

these desired characteristics for known values of the 
area ratio- and Each number are given as figure 2 (pres- 
sure), figure 3 (density), and figure 4 (temperature). 

The mathematical basis for these curves was derived 
as follows: By Bernoulli's equation for adiabatic com- 
pressible flow t 



o 

2 



Pd 



1 P 0 



If - 1 Pd 



(2) 



From 'the requirement of flow continuity, 



whence 



r 3 - V 3 



Po\ S /A 0 i 

Pa 



Also , 



PoV 
7P 0 



and for isentropic flow 
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By substitution of these relations into equation (2) in 
order to eliminate the velocity and density terms, the 
following result is obtained; 




3 M"'* 



7-1 3 



(s) 



for p d /p 0 . in 
terms of M 0 and Aq/A^. The density and temperature 



' This equation' is s oived. in figure 2 Jfa>-t- 0 
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ratios in the duct can "be computed from the pressure-ratio 
values obtained from equation (3). For isentropic flow 



The use of figures 3 and. 4 permits the density and tem- 
perature ratios to be obtained directly for the known val- 
ues of M Q and A Q /A d . 

The use of the charts (figs. 2 f 3, and 4) is not re- 
stricted to the Mach number M Q and the area ratio 

Aq/Aj, as indicated by the subscripts on the figures. 

Tor example, if the area ratio and Mach number were taken 
as A x /A^ and Mj, respectively, the charts would yield 

the ratios Pi/pdt Pi/pd» an(i T i/ T d* other more 

convenient area and corresponding Mach number might be 
similarly used in place of A 0 and M 0 . The charts may 

also be employed to determine the flow from a large area 
to a small area. The Mach number used with the charts 
must always be the Mach number corresponding to the 
smaller area* 

Station 1 . - The entrance area A} should be designed 
to allow the desired mass flow to enter the duct at an 
inlet-velocity ratio that will permit efficient external 
flow as well as efficient internal flow. In the investi- 
gation reported in reference 4, velocity ratios of 0.3 to 
0.6 were found to fulfill these requirements for circular 
openings of suitable shape located at the stagnation point 
of a body. Higher inlet-velocity ratios will make it dif- 
ficult to expand the internal flow efficiently, and lower 
inlet-velocity ratios will generally make the external 
flow unsatisfactory for high-speed airplanes. . If the in-* 
let velocity v x is selected on the basis of these con- 
siderations, the inlet area is given by 




and 
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A, .--JL- (4) 

In most applications a sufficiently accurate value of A a 
will be obtained when p z is taken as equal to p 0 « 

After the entrance area and the ratio & 0 /A>i have 

"been determined, the pressure, the density, and the tem- 
perature at station 1 may be obtained directly from the 
charts if the opening is not located within a boundary 
layer. 

Station 3 »- If the ducting losses between stations 1 
and 2 are a negligible fraction of the free-stream dynamic 
pressure, the static pressure, the density, and the tem- 
perature at station 2 may be obtained directly from the 
charts (figs, 2, 3, and 4) for. the area ratio A 0 /A 2# 

Under most conditions of flight, however, it will be 
necessary' to account for the diffuser losses even in ef- 
ficient systems, The characteristics of the diffuser are 
usually defined by a diffuser efficiency f 



or, preferably,. 



H t * 1 - 



Ah s 



The definition in terms of N the static pressures, as in- 
dicated, can be used only in the reduction of incompress- 
ible-flow data because, at high speeds, the pressure rise 
is affected by factors other than the diffuser shape. 
The definition in terms of the total-pressure loss between 
stations 1 and 2, however, which yields the same value for 
tl^ at low speeds," remains essentially constant with Mach 

number (see reference 5) and should therefore be used for 
all high-speed test data. If the velocity distribution 
is reasonably uniform, the total -pressure loss in the 
diffuser will equal the decrease in static pressuce at 
station 2 belov; the ideal static pressure for zero fric- 
tion loss ; that is, 



Ap 3 



= Ah 2 = ( qi - q 3 )(l - n d ) = qJl - (j^ 




(1 - n d ) 



The assumption of incompressi"ble flow in computing 
the local dynamic-pressure ratio is justified in the 
usual case where the Mach number at the entrance is low 
and where Ah 3 is small. 

The stat ic-pr essure ratio at station 2 is therefore 



P3 /Pa\ <U I^Yl . , N 

P ° \ V °I ideal P ° L ■* *l J 

The ideal pressure ratio is obtained as before from fig- 
ure 2, and the diffuser efficiency can "be estimated 

from test data such as that given in 'references 5 and 6. 
Data on the diffuser efficiency for the model used in the 
present tests are given in reference 5* 

The temperature distribution across the duct at sta- 
tion 2 is not uniform, the temperature being slightly 
higher near the walls where the friction effects predom- 
inate. A satisfactory value is given by 



in which the ideal pressure ratio of figure 2 is used. 
The corresponding value of the density is 




gHT 2 



= 0.000583 




The Mach number of the flow in the duct at station 2 will 
be shown to be an important parameter and is given by 
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M 3 = Mo Z± ±o . Mo Po £o W 



1/2 



Flow changes within the radiator ,- An exact deter- 
mination of the flow changes between stations 2 and 3 must 
include calculations of the flow details within the heat- 
ing device. The following analysis describes the flow 
through a tube or channel of constant cross section in 
which heat is added and in which friction losses occur. 

If no appreciable energy losses are assumed to occur 
at the tube entrances, it is apparent that the flow 
changes . between station 2 and station r 2 (at the tube 
entrance) are determined by the area ratio A r /A 3 and 

the Mach number M r , where A- is the free area within 

the radiator. Figures 2, 3 t and 4 may therefore be used 
to obtain the desired conditions at station r 3 by 
merely changing the nomenclature on these charts. The 
charts require the use of the Mach number 'M rj3 , .which may 

be computed from the known value of ■ M 3 as follows: 




The value of M p ^ is dependent on Pzl? Vz and ^ 2 ^v 2 9 

values which are to be obtained from the charts. A- trial 
and error process must therefore be employed until corres- 
ponding values of M r2t P 3 /pr 3 t and T a/ T r 3 are 0<b ~ 

tained. Two or three trials usually suffice. The pres- 
sure drop of the air entering the tube is given by 



P 3 ~ Pr P s Pa 1 - ~\ <6) 




in which p 2 /p r3 is obtained directly from figure 2, for 
the value of A r;g /A 3 and M r2 as just described. 
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The pressure drop between stations r 2 and r 3 (at 

the tube entrance and exit) is obtained by applying the 
momentum relation between these stations: 

D f . 

^r 3 r 3 *r 3 ^r 3 r 3 • 

2 

where D f is the drag force due to friction within the 
radiator tubes and A r ^ = A r , The pressure drop in terras 
of cu is, thus, 

Apt Pr 8 - Pr 3 D f < Pr 3 v r 3 3 " Pr 3 ^r 3 3 (?) 

q r 3 q r 3 " *r 3 A r 3 " <*r 3 



It is evident from equation (?) that the pressure 
drop in the tube is comprised of two components: ?irst , a 
pressure drop associated with the frictional drag force 
and, second, a pressure drop due to the momentum increase 
in the tube. This second term, frequently neglected, be- 
comes very important when the velocity in the tube is 
large. The drag term is presented in the form of a drag 
coefficient, denoted hereinafter as where 

The basic value of corresponding to the unheated 

low-speed condition Cp can be determined from pressure- 

f i 

drop measurements for this condition: 



a „ A I \Q A J I q „ / 

r^ 2 A r 3 / i \ n r 2 r z f i \ H r 3 / i 



This basic value of the drag coefficient depends on the 
geometry of the radiator tubes and on the Reynolds number. 
A considerable volume of pressure-drop data for radiator 
tubes is available and may be used to compute Cp . 
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When heat is applied, the drag coefficient increases "be- 
cause of the increase in dynamic pressure along the tube* 
If this drag increase is assumed to "be proportional to 
the mean dynamic pressure in the tube, 



0 Df ■ C D, 



or 





(8) 



A 1 + 



Pr. 



In order to evaluate 



CD 



f 



from the known "basic value 
r 3 ' 



the density ratio p^/p,. must "be obtained. 



' f i , " " 3 

This relation is also required in evaluating Ap^/q^ 

(equation (7)). A method for computing p,. /p r exactly 

"a 3 

will now he described. 



. The density ratio will he obtained by solving simul- 
taneously the momentum, energy, and continuity relations 
between the flows at stations r 3 and r 3 . Equation (7) 
showg the momentum relation. The energy equation, relat- 
ing the total energy per unit mass at the two stations, 
is 



v r 3 3 



7 



7 - 1 p. 



m 



2 



7 P *3 



7 - 1 p. 



(9) 



Since the areas at stations r 3 and r 3 are equal. 



p r s v r 3 p r 3 v r 3 



Solving the last three equations simultaneously for 
Pr 3 /Pr 3 gives 
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Pr : 



= 1 + 



H 



c p gmT r . 



+ 0.2M r , 



1 - 



T 3 



+ 7.0 



— 

2 P, 



-1 



(10) 



This relation has "been evaluated for a wide range of values 
of heat input and radiator drag coefficients and the re- 
sults are given in figure 5 for use in design and perform- 
ance work. Values of Pr a /Pr 3 obtained from figure 5 with 
an estimated C Df may "be substituted in equation (8) to 



obtain a more exact value of 



he repeated until the correct values of 
are obtained. 



This process may then 



Pr 2 /Pi 



and 



3" r 3 

Two or three trials are usually sufficient. 



Several interesting conclusions may be drawn from a 
study of equation (10) or figure 5. When the Mach number 
of the flow approaches zero, the value of 



approaches 



1 + 



«p«»*r a 



or 



1 + 



the density 



ratio for heating at constant pressure. As M 



r 3 



become s 



greater, the pressure drop rapidly increases, with the re- 
sult that the density ratio p r /p r _ increases above the 



constant-pressure value 1 + 



H 

c pS ffiT r 3 * 



very great for large values of 
tered with air-cooled engines. 



M 



This effect is 



such as are encoun- 



of 



For each radiator drag coefficient a limiting value 
M r is reached beyond which the solution of equation 

(10) becomes imaginary. This limiting condition occurs 
when the flow velocity at the tube outlets reaches the 
speed of sound, that is. when M- = 1.0. A curve con- 

necting the limiting density ratios is shown in figure 5 # 
When the drag coefficients or the rate of heating is 
high, this critical condition is reached at low values of 
Mj.^ , owing to the large density decrease along the tube. 

No increase in mass flow can occur after this critical 
condition has been reached. It may be inferred from the 
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general knowledge of nozzles that shock waves will occur 
at the tube exits if the pressure drop is increased after 
sonic velocity has been reached* The energy correspond- 
ing to the increased pressure drop is dissipated in these 
shock waves and no increase in mass-flow rate occurs. 
The use of figure 5 as outlined will automatically indi- 
cate when this limiting condition is reached. 

With the density ratio and the drag coefficient eval- 
uated by equations (8) and (10) as just described, the 
pressure drop between the entrance and exit stations of 
the tube is easily calculated from equation (7); that is, 



A Pt r r * 



(11) 



Flow change s across t he radiator .- The pressure 

change between stations 2 and 3 is the sum of the pressure 
drop due to acceleration of the air into the tubes plus 
the pressure drop within the tubes-less the pressure re- 
covery between the tube exits and station 3» The pressure- 
recovery term has yet to be evaluated* If station 3 is 
assumed to be far enough behind the radiator for the ve- 
locity distribution across the duct to be uniform, the 



momentum equation may- be applied between stations 
3 to give 



and 



Pr - P 3 » 2 <U 1 - 



or, as a satisfactory approximation, 



- 1 



(12) 



The over-all pressure drop across the radiator is the sum 
of equations (6), (ll), and (12): 



Ap r =p a -p_ = p 3 1- 



P3 



a -r ; 



CD f + 2 



r r 3 
Pr. 



- 1 



- 2q 




-1 



The pressure drop is usually more conveniently expressed 
as a fraction of the dynamic pressure q 2 , as follows: 
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APr _ 1.43 



Pr 




- 2 



+ 2 



A- 



- 1 



(13) 



For the case of the tubular radiator, simple methods for 
evaluating Pr 3 /p 3 and P TB /pz from figures 2 and 3, 

Pr 3 /Pr 3 from figure 5, and C]) f from equation (8) have 

"been described. Equation (13) for the over-all pressure 
drop is exact except for the approximations made in sim- 
plifying the last term, which gives the pressure recovery 
"between stations r 3 and 3... This. term is a small frac- 
tion of the over-all pressure drop and cannot "be computed 
precisely, in most cases, because the flow distribution 
rarely becomes uniform in the short duct lengths that are 
usually available for the pressure-recovery process. With 
air-cooled engines, for example, the rear pressure station 
is always taken so close to the baffle exits that the 
pressure-recovery term is virtually zero and, hence, only 
the first two terms of equation (13) need be used. An 
example of the application of equation (13) to the case 
of a partly baffled air-cooled cylinder will be given in 
appendix A. 

Simplified method of computing radiator pressure 
drop , - The pressure drop across the heated radiator can 
be determined by a less involved approximate method, ac- 
curate enough for many applications, which does not in- 
volve calculation of the details of the flew within the 
tubes or fins. The method consists of correcting the 
basic pressure drop of the cold radiator to the design 
operating condition by a relation involving the density 
ratio P2/P3 across the radiator. The use of Ps/p 3 
instead of the density ratio within the tubes Pr^/Pr^ 

eliminates the necessity of computing the tube flow and 
is justified by the fact that the two ratios are approx- 
imately equal. The same general procedure followed in 
deriving the exact equation for the pressure drop (equa- 
tion (13)) will be employed in deriving the more simple 
approximate relation. If p 2 is assumed equal to Pr 2 » 



the pressure drop between stations 2 and 
imat ely 



r 3 is approx- 
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Within the tube (cf. equation (?)), the pressure drop is 



Ap t = Ap f + p r v a - p^r, 2 



•3 '3 



*3 1*3 



or approximately 
Ap t = A Pf t 



i 


i + £a\ 


2 





+ . 2 °- a fe) IP 



The pressure rise* iDe.tweeri the tube exits and station 3 is 
(see equation (12)) 



The sum of these three components is the d.esired over-all 
pressure drop; that is, 



Ap r = Ap fi 



I A, 



i- (l + P2 



- 1 



+ 3q; 



r 3/ 



- - 'a (ft - 11 



(14) 



In order to use this relation to compute the effect of the 
density change across the radiator, Apf^ must "be evalu- 
ated from the known (or calculable) value of the over-all 
pressure drop of the c old radiator Ap ri . For the basic 

cold condition p 2 = p 3 , and equation (14) then becomes 



Ap" ri = AP fi + <l3 
whence 



J±3L 
1A 



1-3 



- 2q ; 



A; 



A. 



- 1 
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Ap fi = Ap Pi - q 2 



A 



- 1 



r 2i 



+ 2q 3 



f A 2 



A 



- 1 



(15) 



Thus Apf can easily be computed from the over-all pres- 
sure drop for the cold radiator if. the geometry of the 
radiator is known. The use of equations (14) and (15 ) to 
correct the pressure drop obtained for the cold condition 
(or for a test condit ion with heat) to the design operat- 
ing condition has been found to yield results approaching 
in accuracy the results of equation (13) in cases where 
the Mach number of the flow in the tubes is not large. As 
previously discussed, the last term in both equations (14) 
and (15) should not be included in cases where station 3 
is too close to the tube exits tq permit the velocity dis- 
tribution to become uniform. 



In order to obtain the correct value, of p, 



/p, 



for 



use in equation (14), it is important that compressibility 
effects, as well as heating effects, be. allowed for. A 
simple method of obtaining p 3 /p 3 from figure 5 will now 

be presented. For this method the drag coefficient of the 
radiator will be needed, and this value can be ob- 

r 

tairied from equation (14) as follows (cf. equation (ll)): 



Ap ] 



q 3 A s 



<l3 



- 2 



[Pa 

[P3 



- 1 



(16) 



where Ap r /q 3 is given by equation (14). The flow in the 1 
constant-area duct between stations 2 and 3 is analogous 
to the flow in the tubes and exactly the spe relations 
apply between stations 2 and 3 as have been applied be- 
tween stations r 3 and r 3 . Equations (7), (9), (10), 

and (ll) may thus be used to compute the density ratio 
p 2 /p 3 and the pressure drop across the radiator by merely 
changing the notation to represent the duct flow rather 
than the tube flow. Similarly, figure 5, which is a 
graphical representation of equation (10), may be used to 
obtain P s /p 3 directly. The aforementioned changes in 

notation, applicable to equations (7) to (ll) and to fig- 
ure 5, are as follows: 
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Replace M pg "by M 3 
Replace p r /p r by p a /p. 



u 3 "3 

Replace- ^- — - "by - H 



3 



D f riction . D radiat or A%% P ^ r 
Replace by , or C D by C D 

( lr 3 Ar 3 ^3 A 2 f r 

Figure 5 may be used with the revised notation to obtain 
a. value for P 2 /P 3 ♦ E y use of a trial value of P 2 /P 3 > 

the drag coefficient C-n can be determined from equa- 

r 

tion (16). With this value of Cr) , a new value of 

r 

p 2 /p 3 can be determined from figure 5. This process is 

repeated until corresponding values of Cp^ and P s /p 3 

are obtained. Equation - (14 ) may. be then evaluated for 
the radiator pressure drop. 

For design work, it is frequently necessary to com- 
pute the pressure-drop ratio Ap r /q Q . This quantity can 

be written in the form 




which can be evaluated in terms of the flight Mach number 
and the design area ratio A 0 /A^ in cases where the area 

ratio. is small (say less than 0.2) and where the diffuser 
loss may be neglected. For these conditions, from equa- 
tion (5) or figure 3, 

2.5 




1 + 0.20M ft s 



Therefore , 
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and 




(17) 



Either the approximate value of Ap r /q g given by 

equation (14) or the exact value from equation (13) may 
"be used in equation (17). The' quantities k 0 /k^ and M 0 

are design constants. 

It has "been shown in figure 5 that the density ratio 
across a radiator or air-cooled engine "becomes very large 
as the tube or duct Mach numbers increase, even with low 
rates of heating* Corresponding increases ' in friction 
and in the pressure drop due to momentum increase within 
the tubes have "been indicated in equation (13) or (14). 
These adverse effects due to the density-ratio change 
should he emphasized "because, under present-day operating 
conditions at high altitude, very large Mach mimbers can 
be obtained within radiators and, in particular, within 
the baffles of air-cooled engines. Calculation of the re- 
quired pressure drops will be seriously in error if no 
allowance is made for the effects of density change. 
Equation (13) or (14) or equivalent relations must be 
used. High-speed test data are needed to indicate how 
closely the one-dimensional analysis of this report approx- 
imates the actual flow through an air-cooled engine. 



In order to illustrate the large pre ssure-dr op increase 
that will be required as the operating altitude is increased, 
figure 6 has been prepared for a -typical existing air-cooled 
engine. Figure 6(a) shows the variation of the pressure- 
drop ratio Ap r /q p with altitude for the air-cooled engine. 
Two limiting conditions beyond which the desired cooling 
conditions cannot be maintained- are included. The first 
occurs when the maximum available pressure drop (assumed 
as 1.0 q Q with extended exit flaps) is reached. Flight at. 
an altitude higher than the altitude defined by this limit 
would require the use of a blower. If a blower were used, 
the second limiting condition would occur at a still higher 
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altitude, where sonic velocity would be attained within the 
baffles, A comparison is made in figure 6(h) of the increase 
in pressure dr.op with altitude as predicted by the methods of 
this report with the increase predicted by current methods in 
which no allowance is made for the density change across the 
engine. Details of the calculation and the assumptions 
involved are given in appendix A. 

A study of the effects of increasing the fin area of 
the engine indicates that the cooling difficulties shown 
on figure 6(a) could be. delayed to altitudes higher than 
those indicated if the fin area could be increased. The 
same result could be obtained by the unsatisfactory ex- 
pedient of increasing either the fuel-air ratio Qr the 
head temperature. This fact undoubtedly explains the 
ability of some existing airplanes to operate at very 
high altitudes for limited periods of time. It. is also 
evident that the cooling problem is alleviated by the re- 
duction in power output that occurs if the critical alti- 
tude, of supercharging is exceeded. 

Temperature change across radiator.- The temperature 

rise across the heating device is easily calculated after 
the pressure drop and dens-ity ratio P 2 /p 3 have been ob- 
tained.' The desired temperature-rise ratio At r /T 3 is 
derived as follows: 



!a. - Is. 

*a P 2 P 3 




T a _ A *r 



Pa 




1 




(18) 



where Ap r /q. 3 . is given, by equation . (13) or (14) and 
Po/Pt ia obtained from the chart (fig. 5) using the 'duct- 
flow notation. 
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Station 3 .- If either equation (13) or (14) is used 
to compute the pressure drop and equation (18) is used to 
compute the temperature change, the conditions at station 
3 are given by 



Stat ion 4. - The internal static pressure at station 
4 depends on the shape and location of the outlet. Ef- 
ficient outlets of the type employed in the present tests 
are designed to make the streamlines of both the internal 
and external flows as nearly parallel as possible at the 
outlet; thus, the external outlet pressure "becomes equal 
to the internal pressure. (See references 4 and 7.) The 
value of the outlet pressure for this type of outlet is 
found to be about the same as the pressure existing at 
that station on the body before the opening was added. 
Pressure-distribution data for faired bodies may there- 
fore be used to estimate the outlet pressure for openings 
of the type described. If the walls of the outlet are 
converging rather than parallel, however, the outlet pres- 
sure will be considerably higher than the pressure on the 
corresponding faired body. (See references 4 and 7.) 
With cusped outlets at the tail of a streamline body, as 
assumed for the present analysis, the outlet pressure is 
about 0.11 q 0 above stream static pressure, according to 

reference 5. If an estimate of the outlet pressure is 
made as suggested, the flow conditions at the outlet may 
be computed from 





0.000583 =-a- 
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Y-i 





0,886 



or 



p = 0.000533 jjA 

1 4 



O = p 




1/7 



The size of the opening required to permit passage of the 
design quantity of internal flow can now "be computed from 



A 4 = 

P 4 V 4 



(19) 



where m is known, p 4 has been given previously, and 
v 4 is given- by- 



V 4 = 
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7-1 p r 



1 - ^ 





+ 7.0 



1 - 




0»886' 



(20) 



The assumption of isentropic flow "between stations 3 
and 4 yields accurate results for unobstructed outlet pas- 
sages. In ordinary cases the effects of friction are local- 
ized, in a shallow, "boundary layer near the wall and do not 
affect the main outlet flow. In many practical arrange- 
ments, however, the- outlet passage is complicated Toy turns, 
tiers of guide vanes, and like arrangements, to which spe- 
cial consideration must he given. In such cases, it is 
suggested that the step-by-step process employed in deter- 
mining the flow across the radiator.be used. The drag 
coefficient of the obstruction can usually be estimated, 
and figure 5 may be used to simplify the computations. 
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A less exact "but more simple method of accounting 
for the energy loss due t*o obstructions "between stations 
3 and 4 assumes that the isentropic expansion takes place 
between the pressure p 3 - Ah 3w . 4 and the outlet pressure, 
where Ah 3 _ 4 is the total-pressure loss- occurring "between 

stations 3 and 4; that is, as an approximation, 

p4 ~ P . 3 \ P3 " Ah -V 

A sample calculation using. this method was given in refer- 
ence 8 . v 

For us\ial installations, the maximum duct velocity 
"between stations 3 and 4 will "be attained at-' the outlet. 
It can "be shown that the outlet Mac.h number cannot exceed 
unity. This limiting condition cannot "be reached in ordi- 
nary cases, however, without the aid of an internal blower 
because the outlet Mach .number will not exceed the flight 
'Mach tiumber except in systems where the outlet is located 
in a region of high negative pressure. The addition of 
heat alone. will not appreciably change the outlet Mach 
number because, with the> increase in outlet velocity due 
to the addition of heat, there is a corresponding increase 
in the speed of sound. If the critical velocity were 
reached' in the outlet, the mass-flow rate could still be 
increased by an increase in outlet area. 

Station 5 .- The velocity at station 5' is required for 
the purpose of computing the resultant force due to the 
internal flow. /As the static pressure is equal to p 0 , 
by the general energy equation 

The resultant, force is obtained by substitution of this 
value of v 5 into the momentum relation which t since 

,P 5 = P 0 • is simply 
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D = m(T 0 - v J 



(21) 



A negative result obviously indicates that the resultant- 
force is a thrust. 



The Effects of Heat and Compressibility on the Total Pressure 
at the Outlet and on* the Internal Mass Flow 

Total pressure at the outlet . - In cases where the 

station-to-station method is not employed, the total pres- ' 
sure at the outlet must be measured or calculated in order 
to determine the net drag or thrust due to the internal 
flow. The station-tp-station method of analysis just de~ 

* scribed, although essential in c$lculat ing exactly the . 
details of the internal' flow, does not clearly indicate 
the over-rail effects of heating and compressibility on 
the outiet total pressure. .A more direct, method, which 
is useful in illustrating the effects in question .and 
also in computing the mass-flow changes a'nd the net 
thrust, will now be described. The total-pressure loss. 

. at the outlet is the sum of the loss across the radiator, 
the duct loss, and the pressure rise through the blower, 
if one is present; that is, 




over-all 



(22) 



bl ower 



The total-pressure loss across the radiator is equal to 
the static-pressure drop (equation (17)) less the dynamic- 
pressure rise (for the small values of M 3 usually en- 
countered),; 






From equation (17 ) , therefore, if no. blower is used, 

3 / v 2 • 5 



over-all 




1 + 0.2 v OM 0 3 y 





(23) 



duct 
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This equation is easily evaluated "by use of figure 5 with 
the" duct-flow notation t.o obtain the necessary values of 
p 2 /p 3 . The duct loss should include an estimate of the 

skin-friction lbsses plus the diffuser loss, which, from 
equation (5), is given by 




Tor a constant mass-flow ratio A 0 /A 3f the radiator 
drag increases as a result of the increase in density 
ratio Po/p 3 due to heating. As shown by equation (23), 

the total-pressure loss will correspondingly increase. 
Aside from this well-known \ef feet , however, the addition 
of heat will cause little change in the outlet total pres- 
sure, From the development leading to equation (23), it 
is seen that, for a given mass flow and radiator drag, 
the outlet total pressure will decrease as heat is added 
by the amount q 3 ■- « q 2 . . Owing to the low due t . vel oc i t i e s , 

the quantity (q 3 - ^s^/^o ^e very small- in usual 

cases. This result is of considerable interest in that 
measurement of total pressure in the wake will give no 
direct indication of the propulsive effect of heating. 

When the duct velocities are large, however, and the 
corresponding value of q 3 - q g is appreciable, the fore- 
going development indicates that the effect of heat will 
cause an appreciable loss in outlet total pressure over 
and above the increased radiator drag Iocs, 

For a duct system with moderate internal losses, no 
appreciable drag change of surfaces exposed to the outlet 
fl.ow is to be expected because of the negligibly small 
change in outlet dynamic pressure with heating. The 
Reynolds number of the -outlet flow decreases, "however, 
and some scale effect might therefore be expected, No 
change in the outlet Mach number occurs because the ve-lo- - 
city increase is offset'by an increase of the speed of 
sound in the heated outlet flow, ■ 

Variation of mass, flow ,- The internal mass-flow coef- 
ficient — ™ — or the flow ratio i s a constant for 
Po FV a Po V o 
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a given .system, except for the effects of compressibility 
and heating. (See reference 5.) In order to show these 

effects, the'ratifc P *-V* ' - will be written in terms of M Q 

r o o 



and 



H 



CpgmT 0 



• From the general energy equation, 



2. r y - 1 ?4 



Y-l 

0>4 



- 1 



whence , 



. PoV'V 7 ~ 1 Po *0 



where 



h 4.A Y 



- 1 



• V . 



Ah 
P. 



(24)' 



and 



7M. 



P. 



The density ratio p^/p 0 . in equation (24) can De written 



in terms of the heat added and M Q . 



/y 



If the internal losses are -assumed tc he moderate and the 
pressure, at the radiator to approach stagnation pressure, 
it can be shown fronTthe eriergy 'relatione that '* T 5 is 

related to the heat. added as follows: 



29 



whence 



1 + 



H 

c p gmT. 0 (l + 0.20M 0 2 ) 



£4 
Po 



1 + 



(1 + 0.7PM 2 ) 
± Q 



0.714 



7.74mT Q (l + 0.20M o ) 



(25) 



When the indicated suhst ituti ons for total pressure are 
made, equation (24) "becomes 



/ 








1 7-1 




7 


2 + 7P 4 M 0 2 1 


-1 | 


\ J 









where P 4 /p 0 is given by equation (25), Inspection of 

this result shows that, for fixed values of M 0 and 
Ah/q 0> the mass flow will decrease as heat is applied in 



proportion to p 



The outlet area, therefore, must "be 



increased as he at is added to the flow if a fixed mass- 
flow rate is to he maintained. Equation (26) has "been 
solved for a range of values of Ah/q Q and M Q in order 

to show the effect of Mach number and to facilitate the 
use of this relation in design work. The solution is 
shown on figure 7 for an outlet pressure equal to free- 
stream static pressure P 4 =0. The Mach number effect 

on the flow ratio is seen to he appreciable. Equation 
(26) or figure 7 will yield values of P 4 v 4 as accurate 

as the results of the step-by-st ep process for most ap- 
plications, provided that the internal losses are not ex- 
cessive and that the heat is added at a duct pressure ap- 
proaching stagnation pressure. The values of Ah/q Q and 

p 4 /p 0 required for solution of equation (26) or for use 

with figure 7 are obtained from equations (23) and (25), 
re spect ively , 
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If the flow is considered to be incompressible except 
for the effect of heating, equation (26) can be reduced to 



£±^ = ( P -±) ll - f± - 7 | (27) 



as M 0 approaches zero. This relation was given in ref- 
erence 7, The curve on figure 7 . f or M 0 * 0 is a plot 
of equation (27) for the outlet pressure coefficient 
P 4 85 0, Because of its simplicity and also "because equa- 
tion (26) tends to become inaccurate at very low Mach num- 
bers, equation (27) should be used for M 0 < 0.2 in pref* 

erence to equation (26) in cases where figure 7 is not 
used. Figure 7 applies for only one value of the outlet 
pressure coefficient P 4 - 0. It may be used for values 

of P 4 other than 0 by correcting the re suit obtained 

from the figure to the required value of P 4 according 

to equation (27). 

Direct calculation of the outlet area /- Calculation 
of the size of outlet opening required to pass the design, 
mass flow is one of the most frequently performed opera- 
tions in the design of an internal flow system. If the 
step-by-step process has been carried out-, the outlet 
area may be calculated . from the values of p 4 and v 4 

as already shown (equation (19)), In general, however, 
the computation of A 4 without the detailed intermediate 

steps is desirable. Equation (25), with equation (26) or 

(27), or figure 7 can be used to calculate ? 4 Z 4 , which 

ro o 

permits the outlet area to be obtained directly: 



or 




(28) 
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Effect of an operating propeller.- If the inlet open- 

ing is located in the slipstream, the velocity, the pres- 
sure, and the density in the slipstream should be taken 
as the initial conditions and the same procedure followed 
as for the condition without propeller in arriving at the 
conditions in the duct. In cases where the inlet is im- 
mediately "behind the propeller, it. will he found conven- 
ient to take the flow characteristics in the inlet open- 
ing as the initial conditions rather than the conditions 
ahead of the propeller. The area-ratio charts (figs. 3, 
3 , and 4 ) may he used after substitution of A x and M x 

for A 0 and M 0 . 



The Conversion of Heat Energy to Thrust Power 

The effects of the addition of heat in increasing the 
radiator drag and in reducing the internal mass flow have 
"been discussed. Both of these effects cause changes in 
the drag due to the internal losses, "but neither effect 
results from the recovery of mechanical energy from the 
heat added. In the study of the propulsive effect due to 
heating, therefore, it is desirable to eliminate these un- 
related drag changes. This elimination can "be accomplished 
"by assuming, for comparison with the heated condition, a ■ 
hypothetical condition without heat in which both the 
radiator drag and the mass-flow rate have the same values 
as in the heated condition. This assumption not only 
simplifies the analysis of the recovery process but al so 
greatly facilitates calculation' of the net propulsive 
force, as will presently be shown. In the following anal- 
ysis this assumed unheated condition will be indicated by 
the use of primes with : the symbols. 

The net drag due to the f 1 ow t hrough an int ernal sys- 
tem is equal to the decrease in the momentum of the in- 
ternal flow, provided that the momentum i s measured at 
stations where the static pressures are equal to free- 
stream pressure. 

Prom equation (21 ), 

D = m(V 0 - v 5 ) 

If the system produces a net thrust, this thrust is given 
"by the increase in momentum of the internal flow: 



32 



T = -D s= m(v 5 - *V 0 ) 

The drag decrease (or thrust increase) resulting from the 
recovery of mechanical energy from the addition of heat 
is 

AT = (Thrust with heat) - (Thrust without heat 'but 
with same mass fl-ow and same radiator drag 
as with heat ) 



or 



AT = m(v 5 - v 5 1 ) (39) 



The increase in kinetic energy of the internal flow due 
to the recovery process is thus 

A(K.E.) - £ (v 5 2 - v 5 i 2 ) (30) 

If the efficiency of the process is defined as 

E. iS> recovered 
Total heat added 



Then, from equation (30), 



m(v 5 3 -v 5 > 3 ) 



€ H = 2 JH 



(31) 



The over-all efficiency of conversion of heat energy 
into thrust power is the ratio 

c _ Useful thr u st power _ ATV o _ m ^ v 5 ~ ( 32 ) 
T Heat added JH JH 



This equation can he written in terras of the heat-cycle 
efficiency (equation (31)) as follows: 

= € H c M • (33) 
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The quantity e M represents the efficiency with which the 
kinetic energy recovered from the heat input is converted into 
thrust power. This efficiency is analogous to the ideal efficiency 
of a propeller and, in the same way, accounts for the energy left 
in the wake when the wake velocity is different from V 0 . It will 
"be found convenient in the succeeding analysis to express $n in 
terms of the internal-drag coefficients with and without heat. From 
equations (21) and (33), expressing the drag and mass flow in coeffi- 
cient form, 

e M - — t~t — — w 

Efficiency of the heat -recovery cycle .- If the duct velocity at 
the radiator section is very small in comparison with the flight 
velocity, the static pressure at which heat is added will approach 
stagnation pressure and the pressure change across the radiator will 
approach zero . Under these conditions the heat cycle consists of an 
adiabatic compression from the stream static pressure to stagnation 
pressure, an addition of heat at constant pressure (p 2 a* p r » p Q ), 
and an adiabatic expansion tack to p Q . The continual admission of 
air at pressure p 0 and density p Q is equivalent to cooling the 
outlet flow "back to p 0 at constant pressure p 0 and this step 
completes the cycle. This process corresponds to an ideal heat 
engine cycle, the efficiency of which is, from elementary thermo- 
dynamic relations, 

1- 7 

e _ Kin etic energy recovered - ± - 
H Total heat added 

The assumptions on which this formula is "based are justified in the 
case of the internal flow system only when the duct velocities are 
very small by comparison with the flight speed and when the friction 
losses in the radiator or combustion chamber are negligible . 
Operating conditions exist, however, in which relatively high duct 
velocities occur, and an expression for the efficiency applicable 
to such cases will now he derived. The heat-cycle efficiency will 
he shown to depend not only on the compression ratio hut also on 
the Mach number of the duct flow and on the rate of heating. 
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. From the energy equation (cf. equation (9)), 

5 ., g 5 , - = f ; -,Jgc p (T 5 - T 5 ') ..... 

By equation (31).^ therefore, . 

* H = 1 - -"V (T 5 - V) . (36) 

The values of T 5 and T 5 1 ' will depend on the outlet pressure 
and on the conditions at station. 3,, which in the general case will 
depend on M 2 and on the radiator or combustion-chamber losses. 
The previous analysis has indicated that when, a radiator of con- 
ventional design is present, it. .is impossible to obtain high Mach 
numbers in the duct at station- 2 because the flow will choice within 
the . radiator at the tube exits, .for relatively low values of M 2 . 
Since it is desired in the present analysis to show the effects of 
relatively large increases in on the efficiency of the recovery 

cycle, it will "be . assumed that no . radiator obstructs the' duct • The 
addition of heat will "be assumed to take place by combustion of 
fuel in the duct "between stations 2 and 3 "by means of a "burner 
having negligible friction losses.. This ideal arrangement is being 
approached in the ram-jet types of propulsion system. The effects 
of the increase in mass flow due to addition of fuel are of 
secondary importance and will be neglected. From these assumptions 
T 5 1 can be related to T 2 as follows: 



1-7 




The pressure p 3 is lass than p 2 by the value of the momentum . 
increase due to the . addition of heat; that is, inasmuch as A 2 = A 3 
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Also, from the equation of state , 

Ta = T * p 3 p 2 



arxd, substituting for p 3 /p 




The heat-cycle efficiency is, therefore, from equation (36) 




This equation for c? H is plotted in figure 8 as a function of 
compression ratio for various values of the heat input factor 
and M 2 . The term p 2 /p 3 in equation (37) is a function only of 
these variables and was computed from equation (10) for Cj^ ^ 0, 

using notation appropriate to the duct rather than the radiator 
tube 

The Mach number at the entrance to the combustion chamDer, Mg, 
ie shown in figure 8 to have a large effect on the heat-cycle efficiency. 
When M 2 is very small, approaching zero in value, the heat-cycle 
efficiency approaches the value given by equation (35), regardless of 
the heating rate. As increases, however, the value of e H decreases 
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because of the pressure drop due to acceleration of the heated air 
in the combustion chamber. For a given value of M £ , increasing 

the rate of heating causes decreases in for the same reason. 
The limiting or "choking 11 values of are also shown in figure 8. 

If no heat were added, the highest possible value of Mg would 
be 1.0. As heat is added, the choking value of Mp decreases 

because sonic velocity is reached at the end of the combustion 
chamber (M 3 = 1.0) for values of M 2 less than 1.0. At the 

H 

high rates of heating (for example, at ' gIQ ^r = 5*0, which 

corresponds to a temperature rise of the order of 2^00° F in the 
combustion chamber), the choking value of M 2 becomes very small, 

only 0.196 for the example cited. This effect is extremely important 
because it dictates the minimum possible area expansion A 2 /A 0 for 
a given heating rate. 

The flight Mach numbers required to produce the compression 
ratios of figure 8 (assuming ideal stagnation pressure is obtained 
at station 2) are noted on the figure. At the higher supersonic 
speeds, the ideal thermal efficiencies for low values of M £ 
approach those obtainable in internal -combustion engines. This 
result would be expected in view of the fact that the compression 
ratios are of comparable magnitude. 

Because e H increases with compression ratio, the inlet 

ducting and dif fuser losses should be ac low as possible . If an 
internal blower is not used, high pressures in the heating section 
of the duct can be obtained only when the area of this section is 
greater than the initial area A Q of the internal -flow tube. 

(See fig. 2.) This location of the heating device in a high- 
pressure, low-velocity section of the duct is also ideal from the 
standpoint of minimizing the losses in the heating device. 

The friction loss within the heating device, Ah r , has the 
undesirable effect of reducing the compression ratio. For reasons 
of simplification this loss was neglected in the evaluation of e^, 

equation (37) , figure 8. If the principal part of. this loss occurs 
upstream from the region in which heat is added, it may be accounted 
for by using p^ - Ah r instead of p 2 a3 the effective pressure 

at the entrance to the heating region. In the case in which the 
friction loss is intimately involved in the heating process as in 
a tubular radiator, the actual pressure drop during heating will 
be greater than assumed in equation (37) . As an approximation in 
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this case, it may "be assumed that the effective pressure is p 2 - -5— 

at the entrance to the heating chamber, and this value may then he 
used in equation (37) or figure 8 to obtain In "both of these 

examples, the value of Ah r should include estimates of the heating 
and compressibility effects on the friction loss if low-speed data 
for unheated devices are used in the calculation. 

Calculation of the net internal drag (or thrust) .- j n per- 
formance studies, it is desirable to compute the resultant drag, 
or thrust, due to the heated internal flow without going through 
the step-by-step duct calculations previously discussed. The 
method to be described is exact within the limits of the assumption 
that the outlet flow is uniform. Very few steps and only simple 
slide-rale operations are involved. 

The net force due to the internal flow (equation (21)) is 

D = m(V Q - v 5 ) = m(V 0 - v 5 ') - m(v^ - v 5 <) (38) 

Equation (29) shows the significance of the two terms of eaua- 
tion (38) to be: 

(Net force due to internal flow) = (internal drag 
without heat) - (Drag reduction due to heating) 

Each of these components can be computed without difficulty, and 
the separation of the two is the basis for the present method. 

If the flow were incompressible, the drag corresponding to 
the unheated condition could be very easily computed from 




or 
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where Ah/q 0 is obtained from equation (23) for values of mass- 
flow coefficient and radiator drag corresponding to the heated- 
flow condition. For compressible flow, unfortunately, the drag 
calculation is laborious. A relation exists, however, "between the 
true value of Cj) for compressible flow and the value computed 
from the simple incompressible -flow formula, equation (39) . From 
equations (21) and (39), 




This equation involves only the total -pressure loss Ah/<1 0 
and the velocity ratio v 5 /V Q , which can be defined in terms 
of Ah/q Q and M Q from equations presented in appendix B. A 
plot of this relation is given in figure 9 and may be used to 
obtain the true value of the internal drag for the unheated con- 
dition using the value of the low-speed drag coefficient computed 
from equation (39); that is, 



CD'-cvijga. 




the quantity I 5— 1 being taken from figure 9. 



The drag for the unheated condition having been already 
obtained, only the drag, reduction due to heating (second term of 
equation (38)) remains to be computed. From equations (32) 
&n<l' (33) the drag reduction, or thrust, due to heating is 
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AD = 



or 



AC D = 



(40) 



In summary, the net or resultant force (taken as positive in 
the drag direction) is given by the following formula: 



C D = Cd' - AC D = 2C 




(41) 



in which the parameters involved are computed as follows: 
Ak/'lo equation (23) 

^C D /Ci>J) ' from figure 9 

€ jj ' from figure 8 or equation (35) or (37) 

€ M from equations (34) and (4l) by trial. 

The separation of the propulsive effect of heating from 
the internal drag due to frictional effects is useful in several 
other connections . In the correction of wind-tunnel data obtained 
without heat, for example, the right-hand term of equation (hi) 
may be used to evaluate the drag reduction due to heating. It is 
obvious that model tests should be ma.de with values of Ah/q Q 

and C corresponding to the heated condition. A further applica- 
tion of the separation principle will be made in part II, where 
the thrust force due to the added heat is separated from the 
net force due to the internal flow. The results analyzed on 
this basis, are significant as applied to any internal system. 
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II . HIGH-SPEED WIND-TUNNEL TESTS OF A MODEL RAM— JET 

PROPULSION SYSTEM 
APPARATUS AND METHODS 



The Langley 8-foot high-speed wind tunnel, in which ^ 

these tests were carried out, is a closed-throat, circular— ^ 

section, single-return tunnel with airspeed continuously vm 

controllable for speeds from about 75 to 600 miles per ^ 
hour. The turbulence of the air stream is unusually low 
but is somewhat higher than in free air # 

The model tested consisted of a 15 . 6~inch-diamet er 
body of revolution housing a 1 60-kilowatt heating device 
and mounted on a 9-percent -thi ck , 20-inch-chord , low-drag 
airfoil that spanned the jet (fig. 10). The air inlet 
was at the nose of the body with an expanding duct sec- 
tion back to the radiator. The air outlet was at the 
tail. The general model layout is shown in figure 11 
with the ordinates listed in tables I and II. 

The external shape of the wing-fuselage combination 
was designed to have a critical speed in excess of the 
maximum test speed. In order to meet this criterion, in- 
let and outlet openings designed according to reference 
4 were utilized, and a 9-percent-thick 20-inch-chord wing 
of airfoil section NACA 6S-009 was required. The inlet- 
and outlet-opening shapes of reference 4 had the further 
advantage of permitting wide variations in the rate of 
internal flow and in the inlet- and outlet-velocity ratios 
without causing appreciable ext ernal -dr ag changes. 

Three cusped-type openings at the tail (figs. 10 and 11) 
were employed to vary the internal mass-flow coeffi- 
cient C from about 0.07 for tiie largest siae (tail A) 
to about 0.025 for the smallest rjize (tail C). Tail D 
was similar in size to the intermediate cusped outlet 
(tail B) but had straight, slightly converging walls. 
Tail E was a partial-annular outlet located at the 56- 
inch station on the body and designed according to refer- 
ence 7. Tails B f D, and E had approximately the same 
flow coefficient, namely, 0.052. 

The internal system was designed to house the radi- 
ator and survey equipment with a minimum of duct losses. 
A more detailed description of the model and design 
considerations is given in reference 5. 
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Radiator 

The heat input to the .internal air flow was provided 
"by a 160-kilowatt , 220-volt , three-phase electric resists 
ance radiator designed for a maximum operating temperature 
of 1500° F. The radiator consisted of six elements of 
1V 2 inch X 0.020 Nichrome V furnace strip (radiating sur- 
face! 32 sq ft) woven on a framework of reinforced asbestos 
millboard supports ('figs. 12 and 18). Radiator tempera- 
ture was measured by four chr omel-alumel thermocouples 
welded to the elements. 

Power input to the radiator was controlled by eight 
transformer voltage taps spaced to provide approximately 
uniform power increments. Total power input was measured 
by a standard three-phase . wattmeter. 

Survey Equipment 

Pressure surveys .- Pressures ahead of and behind the 
radiator were measured, respectively, by an 8-tube and an 
11-tube static- and .total-pressure rake spanning the duct. 
(See figs. 11 and 13 for installation details.) All pres- 
sure tubes were connected to a standard manometer and the 
pressure readings were photographically recorded. A 52- 
tube removable survey rake of static- and total -pressure 
tubes was installed at station 4 for measurement' of the 
internal and external pressures at the tail outlet. The 
rake tubes were supported by two 2 V 4 -inch-chord blades 
(NAOA section 16-010) set at right angles to each other 
and extending beyond the outlet walls (fig. 14). The 
tail rake was supported by a lV 8 -inch out side-diameter 
steel tube housing the pressure leads and extending 
through the core of the radiator to the wing duct (fig. 
11). 

Temperature measurements ,- All temperatures were 
measured by calibrated thermocouples connected to a sen- 
sitive potentiometer. Air temperature ahead of the radi- 
ator was measured by two shielded thermocouples located 
in the front thermocouple rake. (See fig. 15 for rake 
details . ) 

Shielding of the junction prevented false air- 
temperature readings due to radiation effects. Air- 
temperature rise across the radiator was measured directly 
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by two thermocouple rakes containing six junctions in 
separate compartments (fig. 16). The thermocouples in 
the front rake were the cold junctions and the correspond- 
ing thermocouples in the rear rake were the hot junctions. 

Temporary apparatus . - During preliminary calibration 
runs, additional apparatus was installed in the model to 
supplement the fixed equipment, A 14-tube total-pressure 
rake spanning the duct was mounted "behind the radiator to 
aid in determining the duct velocity profile. From these 
data the energy added to the air stream was determined "by 
an integration of the mass flow and the air-temperature 
rise. 

The heat lost through the walls of the tail was com- 
puted from measurements of the wall temperatures and local 
skin-friction coefficients. A standard 5-tube static- and 
total-pressure "mouse" was used to survey the boundary 
layer of tail B in conjunction with local wall-temperature 
measurements by thermocouples. 

TESTS 



Drag tests with each outlet were carried through a 
Mach number range of 0,36 to 0,75 (approx, 200 to 560 mph 
in the 8-ft high-speed .tunnel ) . The tail rake (fig. 14) 
was removed during the force tests. Internal pressures 
at stations 2 and 3 were recorded simultaneously with the 
force data, as was the po\*er input to. the radiator. The 
range of power input was from 0 to 160 kilowatts. 

Tests of the model with the tail outlets A, B, C, 
and D were made with the tail rake (fig. 14) in place and 
were similar in range to the force tests. The object of 
these runs was to furnish data for calculations of the 
heat-energy recovery independent of force measurements 
and also to evaluate external-drag changes under various 
heating conditions. A special run was made with tail B 
with an internal resistance plate that produced a total- 
pressure drop of 0.27 q 0 . 

Supplementary runs to permit an estimation of the 
heat lost through the walls included measurements of wall 
temperature and boundary-layer velocity profiles. 

All tests were made at zero angle of attack. 
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RESULTS 



The methods employed to compute the heat-cycle ef- 
ficiency and the duct— flow characteristics from the test 
data are described in appendix B. 

Figure 16 shows the variation of heat-cycle effici- 
ency with compression ratio as determined from the wake- 
survey data. The experimental results are shown for all 
the- tests and therefore cover a wide range of mass-flow 
rates, heat input, Mach number, and ducting losses. The 
heat-cycle efficiency indicated by the force-test data is 
shown on figure 17 for the same range of variables as 
presented for the pressure measurements. 

Comparisons of the flow characteristics measured at 
key stations in the duct with the characteristics com- 
puted by the methods of part I are shown in figure 18. 
Table III shows a comparison of the test values of total- 
pressure loss with the corresponding values computed from 
the theory. Table IV shows a comparison of experimental 
and theoretical mass-flow ratios. Figure 19 compares the 
experimental variation of mass-flow coefficient with out- 
let density ratio with the variations computed from the 
simplified theory (equation (27)) presented in part I. 



DISCUSSION 



Efficiency of heat-recovery process.- The maximum 
Mach numbers of the duct flow at the radiator section did 

not exceed 0.10. Reference to figure 8 or to equation (77) 

will show that, for this low value of M s f the heat-cycle 

efficiency depends only op the compression ratio p s /p Q ; 
that is t 

6 H « 1 - 



This relation is shown for comparison with the experimental 
data obtained from the wake pressure surveys in figure 16. 
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Virtually all of the heat energy theoretically recoverable 
was actually recovered* The data shown represent a very 
wide range of heating conditions and internal mass-flow 
rates. The maximum air-temperature rise in the tests was 
about 300° F and the flow coefficient C varied from 
about 0.025 to 0.070. 

The recovery process was not measurably affected by 
any increase in ducting loss due to heating because of 
the extremely low initial value of the ducting loss for 
the model tested. With tail B, the medium-size outlet,- 
the loss across the radiator was only 0.008 q Q and the 

loss in the converging duct behind the radiator was 
about 0.007 q 0 for a mass-flow coefficient of 0.050, 

The heat loss through the walls between the radiator 
and the tail outlet is estimated as approximately 1 per- 
cent of the heat added. The heat loss determined from 
the wall-temperature and boundary-layer surveys agreed 
with the theoretical calculations of this loss. Allow- 
ance for this small correction was made in the reduction 
of the data, but there was no measurable change in the 
recovery process with so- small a heat loss. 

The force measurements of the net thrust caused by 
heating provided a means of determining a value of €g 

that includes the external losses associated with exhaust- 
ing the hot air of increased velocity. This value of € g 

is shown in figure 17. As noted irom the definition of 
e H (see equations (PI) and (33)), the over-all efficiency 

given by the force tests was reduced by the mechanical or 
Froude efficiency factor- € M in computing ^ The force- 
test data closely approach the ideal value of €g. The 

results are consistent, however, in indicating small losses 
Comparison of figures 16 and 17 leads to the conclusion 
that the losses must be due to slight increases in exter- 
nal drag from the heated exhaust flow. 

Measurements of the external drag with the tail rake 
were made in only the vertical plane (fig. 14), The re- 
sults so obtained showed no consistent external-drag 
changes with heating. In any case, such changes in ex- 
ternal drag were very small; a 1-percent change in the 
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drag of the body would account for the small divergence 
between the ideal and the measured efficiency at the 
highest test speed and heat input. 

The heated flow from the annular outlet had as little 
effect on the external drag as did that from the tail out- 
lets. This result does not seem surprising when it is 
considered fehat the dynamic pressure in the heated flow 
was virtually equal to that of the cold outlet flow. The 
increased velocity is offset by the decreased density so 

that there is little change in l.p 4 v 4 2 and correspondingly 

2 

little change in external drag. (See part I,) 

The efficiencies obtained with the duct pressure re- 
duced by 0.27 q 0 by means of the baffle plate approached 

the predicted efficiency for the measured compression 
ratio as closely as in the runs with low internal losses. 

Figures 16 and .17 indicate that the use of the ideal 
relation for the heat-cycle efficiency is justified for 
design calculations. If the duct-flow Mach number M 2 

is low, say less than '0.10, the heat-cycle efficiency is 
adequately given by the simple relation of equation (35), 

. - 1 ^ 7 

Otherwise, the more exact relation (equation (37)) must 
be employed, in which case eg may be obtained for the 

particular values of U 2 . and L ... involved. 

; c p gmTp 

The data accumulated in the tests are useful in veri- 
fying the relations, developed in part I, concerning the 
calculation of the flow characteristics at key stations in 
the duct , the effects of heat on the over— all flow charac- 
teristics, and the net propulsive force , Comparison of 
the computed and the measured values of these characteris- 
tics will now be made in t.he foregoing order. 
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Comparison of calculated and measured int ernal -f 1 ow . 

characteristics .- The int ernal -flow characteristics at the 
key stations in the duct were computed from the equations 
of part I. on the basis of the following constants, which 
-correspond to one of the test conditions with tail B: 

H, 3tu per second 68.2 

miles per- hour 456 

feet per second 669 

p 0> pounds per square foot ............ 1653 

p 0 , slug per cubic foot 0.001S70 

C D . . 75 

r i . 

^ il / ( lo^ uct » including .diffuse? loss ... .... . . 0.03 

n d 0. 93 

p 4 0.11 

A 4 (outlet area), square foot , .0.0547 

The computed mass-flow coefficient is 0.0399, as" 

compared with- the measured value, 0.0401. The area, . 

A Q = 0.0403 square foot, whicl& corresponds to the computed 

mass-flow coef f ici ent , was used in calculating the pres- 
sure, the temperature, and the velocity at each station in 
the duct according to the equations of part I. The results 
are compared with the measured data in figure 18 and show 
that the flow characteristics of an efficient system, in 
which the velocity distribution is. fairly uniform, can he 
computed with good accuracy*. 

Ef fects of hea t and com pressibility on the"* int ernal, 
t ot al-pr e s sur e lo ss.- Table III gives a comparison of the 
measured variations, of total-pressure loss .for the tests 
with the internal resistance plate with the results com- 
puted by equation (23). Two rates of heating and a range 
of values of Mach number are shown." The reduction of the 
total -pressure' loss with increased heating is due to the 
decrease in internal mass flow accompanying the. addition 
of heat. For a given rate of heating, the total-pressure 
loss decreases as the Mach number M 0 increases in sys- 
tems with a fixed outlet size. This result is caused by 
the appreciable increase in density of the air at the 
radiator section of the duct as the. Mach number increases. 
With a given outlet opening and fixed rate of heating, the 
mass-flow coefficient does not change greatly as the Mach 
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number advances; the increased density at the radiator 
section therefore resulted in decreased velocity ratios 
through the radiator and in correspondingly decreased 
radiator losses. It will he noted from table III that 
equation (23) predicts both the effects of heating and of 
Mach number very accurately. 

Effects of heat and compressibility on the internal 
mass-flow rate .- As shown by equation (36) t the mass flow 
through a system with a given outlet area will be reduced 
in proportion to the square root of the outlet density, 
provided that the internal losses are- negligible. Figure 
19 compares this theoretical variation with the experimen- 
tal results for all of the tests except the runs with the 
added internal resistance. The agreement is excellent. 
For these tests, the decrease in mass-flow rate was inde- 
pendent of Mach number for a Mach number range from 0.26 
to 0.75. 

When an appreciable internal tot al -pressure loss is 
present, as in most practical installations, the change 
in mass-flow rate no longer varies simply as P 4 V* but 

depends on both the internal resistance and on the Mach 
number M 0 . Table IV shows the flow variation obtained 

experimentally with the resistance plate in the duct and, 
for comparison, the flow variations computed by equation 
(26) and also by the approximate equation (27). The 
exact equation (equation (26)) agrees very closely with 
the experimental results, but equation (27) departs ap- 
preciably from the experimental values at the higher Mach 
numbers. It has been found, however, that the approximate, 
equation (27), which does not account for Mach number 
changes, yields sufficiently accurate results at Mach num- 
bers below about 0.30; 

With an appreciable internal loss, the reduction in 
mass-flow rate with heating is less than the reduction 
indicated by the p V2 relation. The reduced total- 
pressure loss accompanying the addition of heat in. a sys- 
tem with a fixed outlet opening tends to increase the in- 
ternal flow and this effect nullifies to some extent the 
throttling effect of the heat. 

Computation of the net force due to the h eated inter- 
nal flow ,- An example of the sample method described in 
part I for computing the net' force will now be given. The 
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same duct characteristics and the same initial conditions 
assumed in the preparation of figure 18 will be used. In 
addition to the flow characteristics previously listed, 

Ah 

the total-uressure loss — = 0.238 and the compression 

P ' ^o 
ratio ~£ = 1.226 are required. These values were 

Po 

computed from the basic data. 

The low-speed drag coefficient without. heat is (from 
equation (39)) 



C D < = 2 X 0.0401 (1 - J 1 - 0.238) = 0.0103 



The drag coefficient without hert "but at the design Mach 
number M 0 = 0.615 is (from fig. 9) 

C D ' = 0.0103 X 0.790 = 0.0081 

From equation (35) (or fig. 8 for M g * 0), the efficiency 
of the heat cycle is 



€ H = 1 - (1.226) 0,286 = 0.0565 

The drag reduction due to heating is given by equation 
(40). The value of € M required in this formula is ex- 
pressed by eauation (34), which cannot be evaluated until 
equation (40) is evaluated. As a trial value, will 
be assumed equal to 1.0. The drag reduction due to heat- 
ing is then, by equation (40) 



» n _ _ JHe H C M - _ 778 X 68.3 X 0. 0565 X 1.0 _ _ 0.0101 
D ' q Q FV 0 • 441 X 1.009 X 569 

The value of the net force is, finally, 



C D = C D « + AC D = 0.0081 - 0.0101 = -0.0020 
If this first approximation is used for the net force, a 
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value for € M may "be computed from equation (34); that 
is, 



1 - 



/ -0.0020 + 0.0081 
i 4 X 0.0401 



= 1.04 



With this value for e M , the new value obtained for the 
thrust due to heating is 

AC D = -0.0101 X 1.04 « -0.0105 

and the net force is 

C D = 0.0081 - 0.0105 a -0.0024 



Ko further trials need be made in this example "because 
€ M is so close to unity. The measured value of the net 

internal drag for this test condition was -0.0022. The 
slight difference between the measured and the computed 
results indicates the small error involved in the assump- 
tion of uniform outlet-velocity distribution. The fore- 
going method is exact except for this assumption. 

Compar is on of subsonic ram-jet propulsion system with 
internal- combust ion engine .- The compression ratio correspond- 
ing to stagnation pressure in the duct at a Mach number 
of 1.0 is smaller than the values employed in internal- 
combustion engines. The ideal heat-cyclo efficiency for the 
duct system considered in this report (ram jet) is, therefore, 
lower at subsonic spepds than for the internal-combustion 
engine. A number of energy losses, however, occur in the use 
of the internal-combustion engine that do not occur in the 
heated duct flow f These losses may be listed as follows: 

1. Mechanical losses in the moving parts of the 
engine 

2. Power required for cooling-air flow 

3. Increase in airplane drag due to the engine 
installation 



4. 
peller 



Friction and compressibility losses 



in the pro- 
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From a consideration of these factors, a value of 24 per- 
cent for the thermal efficiency of a typical engine in- 
stallation was arrived at. This value does not include 
the induced losses in the propeller slipstream or the com- 
pressibility losses. In the case of the duct flow, the 
present results indicate that the ideal thermal efficiency 
can "be closely approached if the ducting losses can "be 
kept low, Furthermore, the results of references 4, 5, 
and 7 show that installation of the duct with suitable 
inlet and outlet openings can he accomplished without in- 
creasing the airplane drag. A comparison of the two sys- 
tems of propulsi on i s shown in figure 20 . The ideal ef- 
ficiency of the duct system becomes of comparable magni- 
tude, to the thermal efficiency of the engine installation 
at Mach numbers approaching 1.0. 

The power recoverable from the cylinder-cooling flow 
("Meredith effect," reference l) in a typical air-cooled 
engine is also shown in f igure 20 ♦ pressure-drop 
results given in figure 6 were used to compute the com- 
press i on r at i o p s /p 0 . The heat dissipated in cylinder 

cooling was ass\imed to be 0.4 brake horsepower. As shown 
in the figure, at a given Mach number, greater recovery 
is possible at sea level than at altitude because the re^ 
quired Ap/q 0 increases rap idly with increase in altitude. 

The energy recovery for this air-cooled engine installa- 
tion thus amounts t o about. 3 percent of . the brake horses- 
power at sea level at M Q = 0.60, but at 30,000 feet 

there would be no recovery at this speed. For lower 
speeds at 30,000 feet the effect of the heat would cause 
a power loss . . 

SUMMARY OF RESULTS 
Compressibility Effects 

The results of both the analytical and experimental 
parts of this investigation are summarized as follows: 

1. Compression of the air in the duct ahead of the 
radiator results in large increases in pressure, tem- 
perature, and density at. the flight Mach numbers now 
commonly attained. 



51 



2. The Mach number of the flow in the radiator is 
an important parameter governing the flow changes across 
the radiator or air-cooled engine. At high internal Mach 
numbers, the density decrease through the radiator due to 
the heat added is greatly augmented by the effects of 
friction and momentum change* 

3. A large increase in the pressure drop, corre- 
sponding to this density decrease, occurs for a given mass 
flow as the Mach number advances. This increased pressure 
drop results partly from increased friction losses and 
partly from the change in momentum of the air flow within 
the radiator. The magnitudes of these compressibility or 
density-change effects increase with radiator drag coef- 
ficient as well as with the Mach number, 

4. In order .to maintain required cooling of air- 
cooled engines as the altitude is increased, an increas- 
ing volume rate of air flow is required because of the 
reduction in atmospheric density. The Mach number of the 
internal flow therefore increases with altitude and the 
compressibility effects consequently become more severe. 
It appears that the use of a blower to produce adequate 
pressure drops at 'high altitudes will be required for 
some existing engines. 

5. The maximum possible mass-flow rate is attained 
when sonic velocity is reached at any station in the in- 
ternal flow system. Increasing the pressure drop after 
this limiting condition is reached will not increase the 
mass-flow rate; the added energy will be dissipated in 
internal shock waves. 

6. The attainment of sonic velocity is the limit- 
ing possible flow condition regardless of the rate of 
heating or the magnitude of the internal-friction losses. 
Increases in either of these factors, however, reduce the 
mass-flow rate obtained at 'the' limiting flow condition. 

7. In conventional cooling-system arrangements, 
sonic velocity .would first be obtained either within the 
radiator at the tube exits or at the. outlet opening. Ex- 
isting air-cooled engines generally employ much higher 
cooling-air velocities than an efficient radiator, and 
the critical flow condition might occur under flight con- 
ditions now attainable. Calculations for a typical air-* 
cooled-engine installation indicated that sonic velocity 
would be attained in the baffle exits at altitudes rang- 
ing from about 40,000 to 45,000 feet if a maximum head 
temperature of 450° F and a fuel-air ratio of 0,10 were 
maintained. 
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Effects of Heating 



1. In usual internal flow systems heat is added in an expanded 
section of the duct where the static pressure is greater than 
atmospheric pressure . Under these conditions the momentum and 

the kinetic energy of the flow in the wake are increased. The 
ratio of the increased kinetic energy, to the heat energy added is 
denoted as the efficiency of the energy-recovery process. The 
momentum increase, which determines the thrust or drag reduction 
due to heating, is related to the kinetic -energy gain by a factor 
depending on the relative magnitudes of the wake velocity and the 
flight speed. 

2. The thermal efficiency of the recovery process depends 
mainly on the ratio of the static pressure in the duct at the 
entrance to the heating chamber to the atmospheric pressure and, 
to a lesser extent, on the Mach number of the flow in the duct, 
and on the rate of heating. Increases in. either the Mach- number 

of the duct flow or the heating rate cause decreases in the efficiency. 

3> An important effect of addition. of heat in a constant-area 
duct is to increase the Mach number of the flow. For heating 
rates corresponding to temperature rises of the order of 2500° F, 
a Mach number of 1.0 is reached at the end of the duct for an inlet 
Mach number cf about 0.2. This critical or choking inlet Mach 
number dictates the area expansion required in the diffuser ahead of 
the combustion chamber. 

1*. The compression ratio, and hence the efficiency, increases 
approximately as the square of the flight Mach number. Maximum 
recovery at a given flight Mach number requires that factors 
tending to reduce the compression ratio be minimized. Poor installa- 
tions having large internal -pressure drops will suffer also from 
decreased gain from the heat-recovery process . 

5 . The power recoverable from the heat added to the cooling 
air (Meredith effect) was found to be about 3 percent of the 
brake horsepower for an existing. air-cooled engine at a flight 
Mach number of 0 .60 at sea level. At 30,000 feet, however, no 
recovery occurred at this flight Mach number, ejid at lower Mach 
numbers the recovery was negative . 

6. The addition of heat in. an internal flow system has a 
throttling effect on the mass -flow rate. The flow for a given outlet 
area decreases approximately in proportion to the square root of the 
outlet density. It is essential that this effect be allowed for in 
computing the required outlet -opening size. 
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7« The outlet dynamic pressure p v 2 does not change 

appreciably with heating because, in most cases , the reduction in 
density is nearly offset by an increase in the outlet velocity 
term v * . The drag of bodies in the wake, therefore, would 

not be changed due to heating, except for possible scale effects 
associated with the reduced Reynolds number in the heated flow. 



Charts and Equations 

For convenience, the formulas and charts showing the principal 
heating and compressibility effects are listed as follows: 

(a) Pressure, density, and temperature ahead of radiator: 
figures 2, 3, and k, and eq.ua tion (5) 

(b) Density decrease within the radiator and across the 
radiator: figure 5 

(c) Pressure drop within the radiator: equation (11) with 
the aid of figure 5 

(d) Pressure drop across radiator: equation (13) or (Ik) 
with the aid of figure 2 arid figure 5 

(e) Effect of density change across radiator on the pressure 
drop: equation (1*0 

(f ) Temperature rise across radiator: equation (18) 

(g) Required outlet area: equation (28) with the mas3-flow 
ratio given by equation (26) or (27) or by figure 7 

(h) Ideal efficiency of heat -recovery cycle: equation (35) 
or (37) or figure 8 

(i) Thrust due to the addition of heat: equation (k-0) 

(j) Net force due to heated internal flow: equation (1*1) with 
the aid of figure 8 



5* 



CONCLUSIONS 



1. The flow details in an efficient internal flow system with 
addition of heat can he computed with good accuracy by the one- 
dimensional compressible flow relations presented in part I . 

2. The actual thermal efficiency obtainable from a ram- Jet 
type of propulsion system is virtually equal to the ideal thermal 
efficiency. The thrust available from ram- jets can thus be calcu- 
lated with good accuracy from simple theory if the rate of heat 
addition to the flow is known accurately. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 21, 19^2. 
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APPENDIX A 

CALCULATION OF PRESSURE CROPS REQUIRED FOR COOLING 
A TYPICAL AIR-COOLED ENGINE AT HIGH ALTITUDES 



Baaic experimental data ,- Test data used to establish the 
"basic sea-level air requirement and pressure drop were obtained on 
a Pratt & Whitney R -2800-27 engine installed in a nacelle in the 
Langley 16-foot high-speed wind tunnel. The operating condition 
and the basic data obtained were as follows: 

Brake horsepower at 2kO0 rpm (normal rating), assumed 

constant with altitude 1600 

Fuel -air ratio 0.10 

Air temperature, T 2 

°F- 100 

°F abs , . , 56O 

Maximum head temperature 

°F ; . . . 4 5 0 

°F abs. ........ # 910 

Average head temperature 

°F U05 

°F abs 865 

Mass flow of cooling air 

slugs/sec ' I.I65 

lb/sec . . . . . 37.5 

Pressure drop Ap p , in. of water , 15 

Frontal area A 2 , sq ft . 17. 

Equivalent leai area ratio A r /A 2 . 0 -li 

Variation of mass flow of cooling a i r with altitude .- The mass 
flow required to maintain a fixed head temperature depends on a 
number of variables, but a discussion of the details of the vari- 
ation with altitude is considered "beyond the scope of this paper. 
A relation currently in use for determining the mass flow will there- 
fore be used without explanation of its derivation. If a super- 
charged engine with inter cooler is assumed, this relation for the 
mass-flow variation is 



m altitude ~ %ea level 



/<T rn\ 1.5 

„ V head 2) sea level 



\ 1'5 



head vT 2j altitude 
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The mass -flow calculations will be carried out for values of the 
flight Mach number 1^ equal to Q.k, 0.6, and 0.8 through the 
altitude range. As pointed out in part I, the temperature ahead 
of the engine T 2 depends on the flight speed as well as on the 
atmospheric temperature. For Army Summer Air, the value of T 2 
wa3 determined from figure k . For M 0 equal to 0 .6, the mass flows 
computed at altitude "by the foregoing relation are as follows: 



Altitude 
(ft) 


To 

(°F abs.) 


T 2 

(°F ab s.) 


Mass flow 
(lb/sec) 


0 


560 


> 600 


46.3 


10,000 




562 


38.O 


20,000 


1*89 


524 


31.7 


30,000 


453 


W6- 


27.0 


1*0,000 


418 




23.4 


50,000 


koo 


U29 


21.9 



The importance of allowing for the temperature rise caused by 
adiabatic compression ahead of the engine is illustrated by the sea- 
level case for which the mass f low was. only 37 «5 pounds per second, 
the compression effect being neglected, as compared with k6 .3 pounds 
per second shown in the table . 



Details of the -pre 3 sure -drop calculations .- The required pres- 
sure drop will be computed as the .drop between stations 2 and r 2 
plus the drop in the baffles. If equation (13) is to be applied 
to an air-cooled cylinder on which the baffling covers only the 
rear half of the cylinder, account mast be taken of the fact that 
approximately half of the heat is added to the cooling air ahead 
of the baffle at virtually constant pressure. Thus, the density of 
the air immediately ahead of the baffle at station 2a is obtained 
from the relation 

e 2 a G pS^ T 2 

The value of H, the heat added to the cooling air, was taken as 
the heat equivalent of 0 .h brake horsepower. Thus, equation (I3) 
for the air-cooled engines becomes 




in which P r2 /p 2 and P 2a /Pr 2 are obtained from figures 2 and 3, 
respectively, for the design values of M r and /A . 

The ratio Pr 2 /Pr^ is obtained directly from figure 5 for the 
design values of M r "and Cj^ . The friction component of the 
engine drag Cj^ is assumed to result only from friction within 
the baffles and was obtained -from the basic engine data as follows: 

(1) The over-all test pressure drop Ap r was 78 lb/sq ft 
and P 2a /p 3 was 1.08. 

(2) Substituting these values in equation (Ik) gives 

78 = Ap fl X J(l + 1.08) + 2 x 1.07 x 51(1.08 - 1) + 1.07(51 - 1) 
or 

= 15-2 lb/sq ft 

V ^-friction q a 1 " 0 ^ 
00 By equation (8), 

As an example, the results of the calculations for a typical 
flight speed and altitude are listed in the following table: 



[Altitude 35,000 ft> M Q = 0 .60] 



Station 

(See. 
fig. 1.) 


P 

(sq f t) 


P 

f..slu<A 
Vcu ft/ 


T 
(°F) 


M 


i 

(sq ft) 


0 


h91 


0.000666 


435 


O.60 


125 


2 


630 


.000788 


466 


.0545 


1.31 


2a 


630 


.000730 


503 


.0567 


r a 


548 


.000637 


^83 


.456 


79.5 


r 

3 


498 


.000543 


514 


.518 



5 8 



The pressure drop across the engine is 

^P r = P 2 - P ro =132 Ib/sq ft 

See figure 6 for pressure drops for other altitudes . 

In order to illustrate the large error involved in the neglect 
of the effects of the density decrease across the engine, the pres- 
sure drop will be computed without regard for these effects. In 
this case, the usual assumption is 



l 2 



= Constant with altitude 



or 



whence 



P 1 

x ^p r X *~g = Constant with altitude 

p 3ea level m 



[ ^altitude \ 2 p sea level 
Pr 9 ^r sea level \m sea leve J P a i ti tude . 



= 9h 1V S( 1 ft 

as compared with the exact value, I32 lb/sq ft. (See fig. 6(b).) 

The approximate method (equation (1*0) for obtaining the pres- 
sure drop across the engine without detailed calculation of the 
flow within the baffles was found to yield approximately the same 
results as the foregoing more exact analysis except when the 
velocity of sound was approached in the baffles. The approximate 
method gives no indication of the attainment of sonic velocity. 
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APPENDIX B 

SEDUCTION OF DATA FOR TESTS 

WITH HEAT ADDITION 
Internal Drag 

No information on the subject of computing; from 
wake-survey data, the drag due to heated internal flow at 
high speeds has "been published. The method to "be de- 
scribed is applicable to the usual experimental .case 
where the rate of heating and the static and total pres- 
sures at the outlet are measured. The momentum relation, 
(equation (21)) is the basis of the calculations. When 
allowance is made for the nonunif ormi ty of the flow, 
equation (21) becomes 



or 



A 4 



It is convenient to write this relation in the form 




(42) 



If the usual assumption is made that the total pressure 
at station 5 (where p K = p„) is the same as the outlet 

total pressure, the various terms in this relation may be 
evaluated as follows: From the isentropic relation, 
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and, also, 





(1 + T>), 



(1 + T>) 



12 



in which the total and static pressures are measured and 
1 + T) is the compressibility factor. A more useful 
formula for 1 +' n than that given in the list of symbols 
involves the total and static pressures as meastired in the 
wake rather than the wake Mach number. This relation is 
as follows: 



1 + r, = l + 0.357 



- h - p 
p(l + ti) 



+ 0.051 



h - p 
p(l + T)) 



+ 0.0018 



\ ? 

h ~ p i 
p(l+ tl) | 




Similarly , since 



h 4 = u s , 




(1 + Tl) 



(1 + Tl) 



The only remaining term required for the evaluation of 
equation (4fi) is (p 5 /P 0 ) • The wake density ratio 



P 5 /p Q depends on. two fact ors: 



(l) the rate at which 



heat is added to the internal flow; and (2) the amount of 
aerodynamic heating, which in turn is governed "by the 



Mr 



and by the internal loss. The term 



Mach number 

(p*/p,J m ay be evaluated by writing the energy equa- 
tion between stations 0 and 5. 
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V* 5 



If the relations p 5 = p 0 , T Q = T 5 (p 5 /p 0 ), p 0 = p 0 gHT 0 . 



H 2 Po v o 3 

C v = J(Y - l) ' and M o = — are substituted in this 



equation, then 



*o P 5 . 



H 



c p gmT 0 



whence , 



P_s N ' 

.Po 




1 » 0.20 M 0 2 (q s /q 0 ) 

l + 0.20 M n 2 + — 2— - 

0 c_gmT, 



.p« 



12 



These relations for (p 4 /P 5 ) VS , (ij^)^, (<l 5 / , 
and (p 5 /p 0 ) 1//S , when substituted into equat ion (4 2 ) .de- 
termine the internal drag corresponding to the point in 
the wake at which the pressures are measured. Integration 
across the outlet opening gives the required value of the 
total drag coefficient due to the internal flow. 



External Drag 

The equations for the internal drag apply for the 
external drag with the following changes: 



(1) 



1 + 0.20 M 0 S (q 5 /q 0 ) 
1 + 0.30 M 0 2 



(2) The integral (equation (43)) is evaluated "between the 
surface of the model and the outer extremity of the wake. 
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Internal Mass-Flew Coefficient 



~ P 0 » Y o " 1 \?6 ?o) 



dA 




dA 



These are the same quantities as already derived in con- 
nection with the internal drag coefficient. 



The heat-cycle efficiency (equation (31)) can he 
very conveniently computed from the difference "between 
the drag coefficient for the cold condition and the drag 
coefficient for the heated condition. From equations (33) 
and (34), 



Inasmuch as the mass flow .is. less for the heated condition 
than for the cold condition, the measured internal drag 
without heat is therefore greater than the desired value 
D T which would he measured if the mass flow did not change* 
For small variations in mass flow, however, it is possible 
to correct accurately the drag for the measured unheated 
condition to the hypothetical unheated condition. The re- 
lation "between the total-pressure loss for the desired un- 
heated condition and for the measured unheated condition is 



Calculation of Heat-Cycle Efficiency 





(Op 1 + op) 
4C 
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The outlet total pressure for the desired condition is, 
then, 




With this value for h 4 , the desired value of C D 1 is 
computed directly from equation (42). This value of C D 
permits computation of €g from equation (44). 

Equation (44) may "be used to evaluate €g either 
for the wake-survey values of Cp and 1 or for the 
force-test values of these drag coefficients. 
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TABLE I.- STREAMLINE-BODY OKU I NATES IN INCHES 



wvx c a ni x i 




went e r 


sect i on 


Streamline tail 


X 


R 


X 


R 


A 


Xh 


0 


0 


0 


6.78 


0 


6.49 


.73 


1.2S 


1.5 


6.80 


3.0 


6.25 


1.46 


1 95 




o • a u 


6.0 


5.93 


2.92 


2.93 


10.5 


6.80 


9.0 


5.52 


4.37 


3.64 


13.5 


6.80 


12.0 


5.06 


5. 83 


4.21 


15.5 


6.73 


15.0 


4.52 


8.75 


5.03 


18.5 


6.71 


18.0 


3.94 


11 . 66 


5.64 


21.5 


6.58 


21.0 


3.35 


14.58 


6.10 


23.0 


6.49 


24.0 


2.75 


17.49 


6.44 






27.0 


2.14 


20.41 


6.67 






30.0 


1.53 


23.00 


6.78 






33.0 


.92 










37.5 


0 



TABLE II— DUCTED-BODY ORDINATES IN INCHES 



Diffuser 


Ducted "body nose ordinates 


Tail 


A 


1 3 


C 


X 


R 


X 


R 


X 


R 


R 


R 


0 

.10 
.25 
.50 
•75 

1.00 

1.75 
7.75 
14.75 
16.75 
18.75 
20.75 

22.75 
25.00 


2.3^ 
2.21 
2.21 
2.21 
2.22 
2.23 
2.28 
2.69 
3.17 
3.34 
3-59 
3.9S 
4.53 
5.59 


0 

.025 

.05 

.10 

.25 
.50 

•75 
1.00 

1.75, 
2.75 
3.75 
4.75 
6.75 
8.75 
10.75 
12.75 
14.75 • 
• 16.75 
1S.75 
19.50 


2.3U 
2.44 
2.50 
2.58 

2.75 
2.98 

3.16 
3.33 
3.73 

4.49 
4.79 
5.29 

5.71 
6.04 
6.30 
6*51 
6.66 . 
. 6.76 
6.78 


0 

3.00 

6.00 
9.00 
12.00 
15.00 

18.00 
21.00 
24.00 
27.00 

28.75 


6.49 
6.25 
5-92 
5.47 
4.89 
4.24 

3.57 
2.95 

2.42 
2.07 

1.985 


6.49 
6.25 

5.91 
5.42 
4.80 

4.09 

3.37 

2.70 

2.17 

1.82 
1.705 


6.49 
6.25 
5. 88 

5.37 
4.71 

3.96 .. 

3. 18 

2.46 

1.87 

1.48 

1.365 


Outlet 
diam. 


3.92 


3.36 


2.68 


Nose radius:' 
0.10 


Outlet 
area 


11.07 


7. 88 


4.64 
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TABLE III 

COMPARISON OP EXPERIMENTAL AND THEORETICAL 
VALUES OP TOTAL-PRESSURE LOSS 
[Tail B, "baffle plate] 



Total-pressure loss, Ah/q 0 



M 


H = 


0 


H * 90 Btu/sec 




Equati on 
. (23) 


Exp eri - 
mental 


Equati on 
(23) 


Experi- 
mental 


0.263 
.347 
.403 
.473 
-.545 
,615 

. .689 


0.299 
. 299 
.298 
.294 
.290 
. 286 
.282 


. 0.301 
.298 
.296 
.294 
.292 
.289 
.236' 










0.218 
. 2 25 
.233 
. 236 
.233 


0.228 
.232 
.235 
. 237 
. .239 



TA3LE IT 

COMPARISON OP EXPERIMENTAL AND THEORETICAL 
VALUES OP MASS-PLOW RATIO 
[Tail B , "baffle plate] 



M 


Mass-flow ratio, P 4 v 4 /p 0 V 0 


H = 0 


H k 90 Btu/sec 


Equati on 
(27) 


Equati on 
(26) . 


Exp eri. - 
• mental 


Equation 
(27) 


Equation 
(26) . 


Experi- 
mental . 


0. 263 
.347 
.403 
.473 
.545. 
,615 
.689 


0.767 
.767 
.767 


0,776 
.786 
.7 90 


0.772 
.781 
.786 












0.672 
.694 
.718 
.736 
.756 


0.621 
. 639 
.657 
.668 
.■67 6 


0. 675 
. 696 
.719 
.734 

' .751 


.767 


. 828 


' .819 



NACA 



Fig.l 




(b) For tubular raj; a for. (c) For baffled air- cooked cy/inden 
Figure 1.- Station designation. 
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Figure 6. - Var/ation of pressure- drop ratio with o/titude for a 
typicoJ air-coo/ed engine. Engine conditions assumed 
are full '-rated power: fuel -air ratio = a/0: maximum 
head temperature = 4 50 °F. 
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(b) Pressure-drop rotio interna of z>ea-/eye/ conc/Zt/ons. Mo= 0.6- 
Figure 6 . - Ccnc/uded- , 
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'(Solution of equation (26) for ^ - O.) 
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Figure 9-" Rat/o of the interncZ-drag coefficient for comprsss/aie f/ow 
to the voJue computed for incompressible fhv/ as a function of 
Mach number. 
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Figure 8. — Vor/at/on of neat-cycie efficiency with 
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various inferno/ Mocn numbers at station 2. 
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/v^t/n? /% - Temperature ratio, internol pressure coefficient, and ye/ocitu ratio 
for toil 3 otM-0. 6/5 . Heat input, 68.2 Btu per sq ft per second 
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F/oure%Q- Comppr/^on of the over-a// efficiency when the arQ//o6/e 
he<?+ is Qcfcfed in a/? interna/- combustion engine and 
/n the cfuct at jfa?nc?t/on pressure. 



